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ABSTRACT 

Relativistic neutron-loaded outflows in gamma-ray bursts are studied at their early 
stages, before deceleration by a surrounding medium. The outflow has four compo¬ 
nents: radiation, electrons, protons and neutrons. The components interact with each 
other and exchange energy as the outflow expands. The presence of neutrons signifi¬ 
cantly changes the outflow evolution. Before neutrons decouple from protons, friction 
between the two components increases their temperatures by many orders of magni¬ 
tude. After the decoupling, the gradual neutron decay inside the outflow has a drag 
effect on the protons and reduces their final Lorentz factor. 


1 INTRODUCTION 

A neutron component in y-ray bursts (GRBs) was proposed 
by Derishev, Kocharovsky & Kocharovsky (1999a,b), and 
detailed calculations of nuclear composition show that free 
neutrons are inevitably present among ejected baryons (Be¬ 
loborodov 2003b; hereafter B03b). Any plausible central en¬ 
gine of GRBs is dense and at least mildly degenerate, which 
leads to its neutronization. During the explosion, the ex¬ 
panding neutron-rich material may undergo nucleosynthe¬ 
sis: neutrons tend to recombine with protons to a-particles 
(Lemoine 2002; Pruet, Guiles & Fuller 2002; B03b). This re¬ 
combination may be successful if the outflow is collimated. 
However, even in an extreme case of a very efficient re¬ 
combination, a significant neutron component is left over 
in neutron-rich outflows because the formation of a-particles 
consumes equal numbers of neutrons and protons. The abun¬ 
dance of leftover neutrons may vary from 10% to more than 
90% depending on the precise parameters of the burst. 

The presence of neutrons changes the theoretical pic¬ 
ture of GRB explosion. Firstly, they may develop a some¬ 
what smaller Lorentz factor than protons. When such a de¬ 
coupling takes place, the last n-p collisions lead to emission 
of observable multi-GeV neutrinos (Derishev et al. 1999a, 
hereafter DKK99a; Bahcall & Meszaros 2000; Meszaros & 
Rees 2000a). Secondly, neutrons decay with time. The decay 
impacts the external blast wave at radii r ~ 10^® — 10^^ cm 
because even an exponentially small number of survived neu¬ 
trons carry an energy much larger than the rest-energy of 
external medium (Beloborodov 2003a). 

In the present paper, we study the dynamics of neutron- 
loaded outflows at early stages of their expansion, r < 
10^® cm, before they are decelerated by an external medium. 
The neutrons decay gradually at all radii r fi 10^^ cm, and 
at small r the decay occurs inside the GRB outflow. The 
decay turns out important at r as small as 10^^ cm. 


Interesting effects also take place at small r < 10^^ cm. 
Neutrons are initially accelerated together with protons be¬ 
cause they are collisionally coupled, and the last n-p colli¬ 
sions before decoupling cause a significant heating. 

We assume in this paper a simple hydrodynamic pic¬ 
ture of expansion driven by thermal pressure and study the 
basic dynamics of a uniform neutron-loaded outflow. We do 
not consider internal shocks or possible dynamical effects of 
magnetic fields. The paper is organized as follows. In sec¬ 
tion 121 we briefly review neutron-free outflows, which have 
been studied previously in detail (see Piran 2004 for a re¬ 
view). In section 121 we derive equations describing neutron- 
loaded outflows and calculate example numerical models. 
Results are discussed in section 0| 


2 NEUTRON-FREE OUTFLOW 

We model the GRB outflow as a steady wind with duration 
toRB, luminosity L, and baryon mass outflow rate M. We 
assume spherical symmetry. This is a good approximation 
also for jets with constant opening angle greater than l/F 
- such a jet behaves as a part of a spherically symmetric 
outflow. The main parameter of the problem is 

p = - 10 ^ - 10 ®. ( 1 ) 

Mc-^ 

We consider three components in this section: radia¬ 
tion, electrons and protons. At small radii r 30Ro, the 
outflow is dominated by pairs, all components maintain 
a common temperature and cool adiabatically. An interest¬ 
ing evolution begins at r > 50Ro where can be neglected. 
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2.1 Opaque stage 

2.1.1 Outflow acceleration 


As long as the outflow is optically thick, electrons, protons 
and radiation behave as a single relativistic fluid. The fluid 
has four-velocity U°‘ = (FpC, Fp/JpC, 0, 0) in spherical coor¬ 
dinates {t,r,9,(j>). The trapped radiation is isotropic in the 
fluid frame and described by the blackbody law with a tem¬ 
perature Tr. 

The total stress energy tensor of the fluid is then given 
by (e.g. Misner, Thorne & Wheeler 1973) 

n-a 0 _ („ rri , „ J2\ U°‘U^ ^ 1 

^ ^ ^ Tr T Pp c J ^2 3 9 5 (2) 

where a = 7.56 x erg cm“® K~^ is the radiation con¬ 

stant and g°‘^ = diag(—1,1, r^,r^ sin^ 9) is Minkowski met¬ 
ric in coordinates {t, r, 9, (f). The plasma contribution to the 
energy density is taken equal to pp (? where pp is the proper 
mass density of protons; small contributions from e and p 
thermal motions and the electron rest mass are neglected. 
The four-vector of baryon mass flux is given by, 

F“=ppt/T (3) 


The outflow dynamics is governed by the conservation 


laws. 


1 d(r^ T*") 
dt r^ dr ’ 

(4) 

dF' 1 dR'^F'') 

(5) 

dt dr 


As long as r < tcRB cFp (which we assume to be satisfied 
thereafter) the outflow may be described as a steady-state 
wind, so that d/dt = 0.^ Then the conservation laws give 


L — Anr^cPpTp 

^^oT/ -|- ppT^ = const. 

(6) 

M = 47rr^c/3pFp 

Pp = const. 

(7) 


The ratio of these two constants is the parameter p (eq. 0. 

The description of fluid dynamics will be complete once 
we specify the evolution of radiation temperature Tr with 
radius. The photon-to-baryon ratio in GRB outflows is ~ 
10® (see e.g. eq. 68 in B03b) and their internal energy is 
strongly dominated by radiation. Energy exchange between 
radiation and plasma practically does not affect Tr and it 
follows the adiabatic law. 



where n = pp/rrip is the proton number density; To and no 
are constants defined at the base of the outflow r = 7?o ~ 
10® - 10'^ cm. 


To = 


3T 


IOttFq/So Rgca 


1/4 


1/4 


IQlO^ 


R, 


1/2 


K, 


( 10 ) 


no = 


; 6 X 10^ 


; Ls2 /■ 


V300/ 


-( 11 ) 


dir Tof3oRQmp c® rj 

where Fq = Fp(i?o) and we assume 

ToPo = 1 . ( 12 ) 

Equations 0 , 0 and © give a closed description of 
the outflow dynamics at the opaque stage of its expansion. 
Combining the equations we And 

Fp(l -1-*) = rj, 

where 


4 aTf 


3 nrUpC^ 


= Xo 


(- 

\no 


\ 1/3 

) 


X0 = ^ -1, 

i 0 


(13) 

(14) 


and obtain the algebraic equation for Fp(r), 

(^-.) (d-A 

(see also DKK99a). 


(15) 


2.1.2 Thermal balance 


The evolution of electron temperature Te and proton tem¬ 
perature Tp in general depends on their energy exchange 
with each other and radiation. It turns out that all com¬ 
ponents of the neutron-free outflow maintain the common 
temperature Te « Tp « Tr during the opaque stage. We 
shall verify this with an accurate thermodynamic calcula¬ 
tion. 

The electron and proton components obey the first law 
of thermodynamics, 

Ht) 

dUi = dQi + {Ui + Pi) — , i = e,p, (16) 

n 

where Ui is internal energy density of component i, dQi is 
heat received by component i per unit volume and Pi is its 
pressure; all the quantities are measured in the fluid frame. 
Ui and Pi are related to temperature T, 

nkTi 


Ui = 


Pi = nkTi, i = e,p. 


(17) 


7i - 1’ 

where pi is the adiabatic index of component i; it equals 5/3 
for nonrelativistic electrons and protons. Then from equa¬ 
tion we obtain the equation for T, 

3 dTi Ti dn 1 dQi 1 


Ti dn 1 dQi 
2 dr n dr kn dt' c Fp Pp 


i = e,p. 


(18) 


where dt' = dt/Vp is time measured in the fluid frame. Elec¬ 
trons exchange energy with photons via Compton scattering 
and with protons via Coulomb collisions, 


dQ, 


, = -nk 

dt' 2 


Tc-Te 


+ 


T —T 

n p. 


where 


TT rUp 


2 meOrclnArie Knif.c'^ nipC 


kn 


kTp 


3/2 


(19) 


( 20 ) 


^ At radii r < tQng cFp the leading and trailing parts of the 
outflow are causially disconnected. Therefore the outflow behaves 
as part of a steady wind despite the fact that geometrically it is 
a thin shell already at r > toits c. 


is the Coulomb timescale (with Coulomb logarithm InA ~ 
15; see Stepney 1983), 


TC 


3 me c 

8 UradUT ’ 


( 21 ) 



















Neutron-loaded outflows in gamma-ray bursts 3 


is the Compton timescale, 


Urad = aT^ 


is the radiation energy density, and Tc is the Compton tem¬ 
perature of radiation. For blackbody radiation Tc = Tr- 
The thermal balance for protons reads 


dQp _ 3^^(rp — Te) 

dt' 2 Tep 

We have neglected here the 
tons and radiation; it is ^ 
radiation energy exchange. 


( 22 ) 

energy exchange between pro- 
(mefmpfl times the electron- 


2.2 Transparent stage 


Electron scattering dominates opacity of the outflow, and 
its optical depth is 


Ue fJT r 



L (7t 

47r rrUpC? Fp rj ’ 


(23) 


where rie = n is the electron number density. At radii r > Rr 
where rr < 1, the outflow is transparent and equation 0 is 
no longer valid. The photon luminosity at this point is (see 
eq.0. 


At = ^RlcPpVlaTfl (24) 

where all quantities are taken at r = R-r- Approximately 
this luminosity escapes to a distant observer as a blackbody 
radiation with observed temperature YpTrlRr)- 


2.2.1 Outflow acceleration 


After the transparency radius, the outflow may still be ac¬ 
celerated by radiation. 

A freely propagating photon at an angle 9 with respect 
to radius satisfies the relation r sin 9 = const and becomes 
more beamed at larger r. A typical photon is emitted at 
angle n/2 in the plasma frame and has initial beaming angle 
9rad ~ l/Ap at r = i?T. At larger radii the beaming angle 
decreases as 


9 


rad 


Rr 1 
r Fp{Rr) 


(25) 


One can define a frame where the freely streaming radiation 
remains approximately isotropic. The velocity of this frame 
is (3rad = cos 9rad uud its Loreutz factor is 


F rad — 


sin 9rad R- 


-FplRr 


(26) 


The streaming radiation will tend to accelerate the plasma 
outflow since Frad > Fp. The radiation flux Fly = L.y/47rr^ 
exerts a force on an electron moving with velocity /3p (see 
Beloborodov 2002, eq. A6), 


^ ^ ( i-Pp \ ( _ Fy\ 

dt c \1 + Pp)\ Ft^J 


(27) 


and the plasma Lorentz factor grows, 

p CT 'r L'y 

dr 16-Kr^mpC^Fp 



(28) 


An easy estimate of the acceleration effect may be obtained 

p4 

neglecting the factor in parenthesis. Then equation lEHJ 

rad 

gives 


r^(r) = FflRr) + 


ScTL-y 

Idnmp 




(29) 


The acceleration at the transparent stage is significant if at 
r ^ Rr the second term on right-hand side exceeds the first 
term. This condition is equivalent to aTr > nnipC? at r = 
Rr , which requires most of the outflow energy be in radiation 
at the moment of transparency. The same condition may be 
expressed in terms of rj (e.g. Meszaros & Rees 2000b), 


Tj > rjrad 


Lar \ 

Fti RonripC? J 


10" Llfl R, 


. 1/4 
0,7 ■ 


(30) 


2.2.2 Thermal balance 


The thermal balance of the outflow at the transparent stage 
is still given by equations ira . cnj and I22II . The Coulomb 
timescale is given by the same equation and the Comp¬ 
ton timescale is given by equation lHU with 


Urad{r) 


Fnr^cFp 


(31) 


being the radiation density measured in the plasma frame. 

We note that the outflow is subject to significant Comp¬ 
ton cooling even at r > Rr because the scattering rate per 
electron is still very high at Rr (it is n-y/n ^ 10® higher than 
the scattering rate per photon). The Compton temperature 
Tc appearing in equation cu represents the effective tem¬ 
perature of radiation observed from the plasma frame and 
is proportional to the average photon energy in this frame. 
Tc may be evaluated as follows. 

The typical photon at r > Rr has angle 9{r) = 9rad(r) 
and constant energy e(r) = e {Rr)Fp(Rr) in the lab frame. 
Its energy in the plasma frame is 


e (r) — cTp)! flpCOs9rad') — ^Fp(l flpflrad) 

The Compton temperature changes with radius as 

Tc{r) ^ e'{r) _ 1 Fp{Rr) 

Tc{Rr) e'{Rr) 2 rp(r) 

where Tc{Rr) = Tr(Rr). Equation 1821 completes the de¬ 
scription of Compton energy exchange at the transparent 
stage. 


1-k 


Kir) 


^Ldir) 


(32) 


2.3 Numerical models 

Numerical models of outflows with rj = 200 and 800 are 
shown in Figures Q and |2] Both models have rj < rjrad and 
the outflow Lorentz factor saturates before transparency. 

Upper panels in Figures 0 and show the evolution of 
the Lorentz factor. As long as enthalpy exceeds rest-mass 
energy (a; > 1), the outflow accelerates with Fp oc r. When 
most of the internal energy is converted to kinetic energy 
{x < 1) Fp tends to a constant asymptotic value Fp/ = rj. 
The characteristic saturation radius Rs is where a: = 1; Fp = 
rjf 2 a.t this radius and 
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Figure 1. Lorentz factor (upper panel) and temperature (lower 
panel) as functions of radius for a neutron-free outflow with 
L = 10^^ erg/s, 77 = 200 and Rq = 10^ cm. Electrons never 
decouple thermally from protons and their common temperature 
is shown by the solid curve. Dashed-dotted curve shows the radi¬ 
ation temperature. 



RADIUS [cm] 

Figure 2. Same as Fig. [Tl but for rj = 800. Tp and Te are shown 
by solid and dashed curves, respectively. The thermal decoupling 
Te > Tp takes place at radii r ~ 10 ^^ — 10 ^® cm. 



Figure 3. Same as Fig. m but for 77 = 8 X 10^. The outflow 
becomes transparent before the acceleration stage ends. 


Rs = 2 -'^^SRo- 


(33) 


Solutions for temperatures Tr, Te and Tp are shown in 
the lower panels of Figures 0 and At very small radii 
where kT ^ me(? = 511 keV, the number density of 
pairs is comparable to that of photons, the energy density is 
(ll/A)aT^ and the temperature is reduced by the factor of 
(4/11)^^'* ~ 0.78; this correction is neglected in the figures. 
As long as the outflow is optically thick, all three compo¬ 
nents maintain a common temperature T via Coulomb col¬ 
lisions and Compton scattering, and T decreases adiabati- 
cally with index 7 = |. The outflow becomes transparent at 
radius 


R-r = 


Lax 


A-jvmpC^r^j ' 


rj < rjrad- 


(34) 


At the transparent stage the plasma is still tracking the 
temperature of (freely streaming) photons: Te ~ Tp ~ Tr ~ 
const until the electrons decouple either from radiation (e -7 
decoupling) or from the protons {e-p decoupling). 

The e -7 decoupling occurs at a radius Re^ where the 
Compton time-scale ('ea. l21ll exceeds the expansion timescale 
'Cexp — R/^p e, 


2 LjOt 2o'tL / 7?o \ 

Stv meC^Tpj: ‘i’KnieC^ Vi^r/ 


(35) 


The e-p decoupling may happen or not depending mostly 
on the value of p. The two possible cases are illustrated in 
Figures Qand |21 

(i) Electrons and protons are still coupled at Re-f {jep < 
Texp) and begin a common adiabatic cooling, Tp = Te oc 

. They will not decouple later because Tep = const while 
Texp keeps increasing. This regime takes place 650 and 

is illustrated in Figure Q 

(ii) Electrons decouple from protons before Re-y - The radius 
of e-p decoupling is found from condition Tep = Texp- Using 






























Neutron-loaded outflows in gamma-ray bursts 5 


Te ~ TriRr) = Tr{RB)(RT/Rs) in equation 1201 for Te; 
(and neglecting the small term kT-p/m-pC?), we get 


R. 


~ (687r mpc"') 


n^T- 


3/2 


(f)' 


(36) 


where Ts = ^TqIt] is the temperature at the saturation ra¬ 
dius. At i?ep < r < Rsj, Te ~ Tr > Tp. At r > Re-, the elec¬ 
trons cool down adiabatically and Tep ~ const while Texp 
keeps increasing. Therefore, electrons and protons eventu¬ 
ally regain the thermal coupling. This regime takes place at 
high r; ^ 650 and is illustrated in Figure]^ 

Figure 13 shows a model with 77 = 8 x 10® > rjrad- In 
this case, the outflow acceleration continues after the trans¬ 
parency radius 


Rt — Ro 


f f^rad j 

\~ ) 


1/3 


1] > Wad, 


(37) 


and Fp saturates at 

/ \ 1/9 

Fp/ « Wad ( ) < V, (38) 

y 'fjrad J 

after a few Rr- Since density is lower than in the previous 
examples, the e-p decoupling happens much earlier than the 
e -7 decoupling {Rep R-,). 

In summary, the thermal evolution of a neutron-free 
outflow is dominated by adiabatic cooling. The opaque out¬ 
flow has a common temperature Tp = Te = Tr which de¬ 
creases as r~^ during the acceleration stage and as 
during the coasting stage. Then, after an intermediate stage 
where the details of coupling between e, p, and radiation are 
important, the outflow is again described by a simple adi¬ 
abatic law Tp « Te oc /pjjg presence of neutrons will 

change this picture. 


3 NEUTRON-LOADED OUTFLOW 


We now consider an outflow with a neutron component and 
denote its initial neutron richness (neutron-to-proton ratio) 
by ^ 0 . All other assumptions are the same as in section|21 In 
particular, we consider spherical expansion driven by radia¬ 
tion pressure. Neutrons and protons are injected at r = Ro 
with an initial density ratio f = ^ 0 , and then ^ evolves with 
radius because neutrons continuously /3-decay into protons. 


n ^ p + e + u. 


The mean life-time of neutrons in their rest frame is rp « 
900 s, and the corresponding mean radius of /3-decay is 

Rp = J cPaVndT ~ 8 X 10 ^® 

where r„ is the neutron Lorentz factor and f3n ~ 1; 
r„/ is the final value of F^ achieved at r ^ ( 10 ® — 
10®)Ro ^ Rp- The neutron population is gradually depleted 
as exp(—r/R/ 3 ) and the n/p ratio (measured in the fixed lab 
frame) evolves with radius as 


r„rin ^ 

TpUp 2 ^ _j_ 


(40) 


where n„ and rip are proper densities of the neutron and 
proton components. 


3.1 Outflow acceleration and n-p coupling 


The total luminosity of the outflow (cf. eq. |3 now includes 
the contribution from neutrons, 


L = 47rr®c 



+ PpC 


+PnTipaC 


const, (41) 


where pp = npnip and pn = nrini„ \ the thermal energy of the 
plasma and neutrons has been neglected compared to their 
rest-mass energy. We allow here the neutron component to 
have a different Lorentz factor F„, which will be close to Fp 
as long as the n-p coupling is efficient. 

The baryon outflow rate is given by 

M = 4 TT r® c (/3pFpPp -I- /3„F„p„) = const. (42) 


The outflow expansion is accompanied by adiabatic 
cooling which determines Tr{r). Expansion of volume can 
be described by the decrease of number density of original 
protons injected at the base of the outflow. We denote this 
density by n and distinguish it from the total proton density 
Up that includes the decayed neutrons. They are related by 


rip 


= n 


1 +Co 

1 +c ■ 


(43) 


The radiation temperature at the optically thick stage obeys 
equation (|HJ- We neglect destruction of neutrons at this 
stage (see section 13.611 and assume in this section ^ = ^0 
and n = np. The difference between rip and n will become 
significant at larger radii comparable to Rp. 

Acceleration of the optically thick outflow is described 
by equations lUTt . iBU and 0 , from which we derive 


dFp _ Fp 2a; dF„ 3^0 

dr r 2a; -I- 3 dr 2x + 2,^ 


(44) 


where 

_ 4 aT^ ( 

* 3 UpinpC^ Vrio/ 


a;o 


>7(1 + Co) 

Fo 


1, (45) 



The first term on the right-hand side of equation 14411 de¬ 
scribes the acceleration by radiation pressure; the second 
term describes the deceleration caused by transfer of mo¬ 
mentum to neutrons. 

As long as the collisional n-p coupling is strong, r„ « 
Fp, equation lHHl yields 


dFp Fp 2 a; 

dr r 2x + 3(1 -I- ^ 0 ) ’ 


(47) 


Acceleration of coupled n and p begins to saturate when 
a; « 1 -I- Co at a radius 


fj 

Rsh — pRo, ^p{Rsb') — 


(48) 


At that point p(,c® ~ aT^ where Pb ^ Pp + Pn- In the case of 
^0 1 , acceleration saturates when x is still large because 

the protons are coupled to a large number of neutrons and 
have a lot of effective inertia. Decoupling from neutrons at 
r ^ Rab would allow the protons to accelerate further and 
reach Lorentz factors Fp ~ p(l -|- ^ 0 ) if the outflow remains 
optically thick (Fuller at al. 2000). 

Description of n-p decoupling will require one more 
equation that specifies momentum exchange between neu¬ 
tron and proton components. Neutrons accelerate because 
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they collide with the accelerating protons. In its rest frame, a 
neutron experiences Trein^an-pC collisions per second^ where 
(Trip « 3 X 10~^®cm^ and 

r.e;=r,r„(l-/?p/3.)«i + (49) 

is the Lorentz factor of the neutron component relative to 
the plasma component. Assuming isotropic scattering in the 
center-of-momentum frame, the mean momentum gained 
by the neutron per collision equals Pp = fiVrei^rei where 
jj, = mpm-n/im-n + nip) « mp/2 is the reduced mass. The 
momentum gained by a neutron during time dt is 

dpn — ^ UpT^^^lUfipPreinipC dt. (^d) 

The corresponding change of the neutron Lorentz factor in 
the lab frame is found from the Lorentz transformation of 
4-momentum cip“ = (0,dp„), 

dr„ = r„/3„-^, (51) 

rrin C 

and one gets 

^ ^ ripFlslPreinnp. (52) 

ar 2 

Equation closes the set of dynamic equations describing 
the outflow acceleration at the optically thick stage. 

The collisional coupling between the proton and neu¬ 
tron components is friction, which inevitably dissipates en¬ 
ergy and heats both components. This heating will be in¬ 
cluded in the plasma thermal balance below (section 3.5). 


3.2 n-p decoupling and transparency 


As long as the n-p coupling is efficient, Tp « r„, the relative 
velocity /3rei 1 may be evaluated using equations iH7l and 


2 rp _ 2x 

2a;-b 3(1-h Co) 

2rp_1_ 

UpOnpr 1 -b {9/8)pbc'^/aTr ’ 


(53) 


The decoupling of r„ from Tp may happen at r ^ Rsb where 
aT^ > pb(?, and according to equation I53> 


or 

/Irei « - r < Rsb- (54) 

'fT^pO'np'R 

The decoupling condition is expressed by setting /3rei = 1 in 
equation ilKHi . The /3rei approaches unity at r .C, Rab if the 
outflow has a sufficiently high rj, 


ri>ri. 


Lg np 

4.8 X 10^ 

(LbA 

47ri?o»TipC®(l -b Co) 

(l+Co)!/" 

{R 0 . 7 J 


.(55) 


Then the decoupling radius is given by 

1/3 

Rnp « „ „ - ' Ro, 


8TTRorjnipC^{l + Co) 


(56) 


^ To the first approximation the rate of n-p collisions < av > 
does not depend on the relative velocity v of colliding particles 
and one may take < av >=< anpC >. 


and the neutron Lorentz factor at decoupling is 


Rnp ,, 7?^" 

Ro ~ 171/3 ’ 


1 ? > V*. 


( 57 ) 


No significant decoupling Tp > r„ happens for 77 < 17 ,. 
This is so despite the fact that 17 < 17 * does not exclude 
PpR < aT^ < pbC? at r > Rab, so there may be enough 
energy in radiation to accelerate protons to Tp > r„/. A 
detailed analysis and numerical models show that if /3rei is 
small at r Rab, it remains small at r > Rab- Therefore, 
to a good approximation, 17 > 17 * may be taken as a true 
condition for decoupling of Tp from r„. Hereafter we focus 
on 77 > 77 *. 

The decoupled neutrons do not affect the remaining e-p- 
7 outflow until a much larger radius where /3-decay becomes 
important (see below). Therefore saturation of proton ac¬ 
celeration and transition to transparency may be described 
as if there were no neutrons. This ‘neutron-free’ outflow has 
luminosity 


L = L - -4^ 

1 -b/ 

and mass outflow rate 


M = 


M 

1+Co' 


(58) 


(59) 


Equations of section 2 then apply if one replaces L L, 
M ^ M and p ^ r) = L/Me?. 

The outflow is still opaque to radiation at the n-p de¬ 
coupling radius. This follows from the small ratio of cross 
sections anplnr ~ 1/20. If cascade initiated by is ne¬ 
glected (see DKK99a and section 3.7 below), transparency 
comes soon after R„p. In this case, one can show that an 
outflow with Tp > r„ (77 > ?7*) becomes transparent be¬ 
fore saturation of Tp. This may be seen from the following 
relation, 


V 

Prad 


1 + Co 1 



where (cf. eg. 13011 


Lar \ 

4:TvRompC^ J 


(60) 


( 61 ) 


?7 > 77* implies 77 > rjrad and hence Tp does not reach the 
maximum possible value Tp^max = 77. The outflow becomes 
transparent and Tp saturates at (cf. eg. 1381 



(62) 


It is not much larger than r„/. The cascade initiated by 
7 r° decay may prolong the opaque stage and increase Tp/. 


3.3 Deceleration by /3-decay 

Next, we consider larger radii r > Rn when the outflow is 
composed of decoupled radiation with luminosity de¬ 
coupled neutrons and plasma. The neutrons are no longer a 
fluid; they retain a mildly-relativistic velocity dispersion ac¬ 
quired at their last collisions at r ~ Rnp, which implies that 
r„ varies by a factor < 2. We will neglect this dispersion 
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and assume that all neutrons have equal Lorentz factors r„ 
after decoupling. 

The neutrons continuously /3-decay in the outflow and 
create a source of protons and electrons moving with respect 
to the plasma frame with a Lorentz factor Frei- This beam 
shares momentum with the plasma, heats it and reduces the 
plasma Lorentz factor Tp. 

The outflow luminosity at the transparent stage may be 
approximated as 

L = 47rr^c [PpVllppC^ hp) + /3„r^p„c^] +L^ = const.(dS) 

We include here the enthalpy of proton component hp = 
Up Pp for completeness. As we shall see the plasma is 
heated to a high temperature at large radii r ~ Rp, however, 
hp remains smaller than PpC^. The enthalpy is related to 
proton temperature by 


hp 


7 

7 - 1 


npkTp , 


(64) 


where 7 is the adiabatic index of the proton component, 
which is between 4/3 and 5/3. From equation I42II . 16311 and 
using the hrst law of thermodynamics leas. 1161 and IT71 we 
derive. 


dFp 

dr 


1 

l-\-Xp{2 - 7 ) 


^(r„-rp) , 2 ( 7 - 1 )^ , 
Rp + r 


-b 


rp 7 dqp FpXp dy 
ppC^ dr 7 dr 


otL^ 

IGnr^mpC^Fp 



(65) 


where Xp = hp/ppC?. The last term describes the radiative 
acceleration (see section 2 . 2 . 1 ). 


3.4 Summary of the dynamical model 


Here dQi/dt' (i = e,p) are the heating rates of electrons and 
protons. 

The thermal balance of protons is signihcantly changed 
compared to the neutron-free case because of two effects: (i) 
frictional heating due to n-p collisions and (ii) heating due 
to /3-decay. 

The rate of n-p collisions per unit volume in the plasma 
frame is 


7lev — TlryTlrjO'- 




where n'„ = FreiUn is neutron density measured in the 
plasma frame. The mean energy dissipated per collision is 
{l/2){Frei — F)mp(? assuming that the relative bulk velocity 
is isotropized in collisions. This gives the frictional heating 
rate, 

^ Trel/Trel “ l)mpC®CTnpnnnp. (67) 


Heating also results from the gradual /3-decay of the 
neutron component. The decay rate per unit volume is 
Lorentz invariant and in the plasma frame (where the decay 
time is Freixp) may be written as 


hp 


n 


n 


r relXp 


n„ 

T0 ' 


( 68 ) 


The decay products form an e-p beam with velocity Vrei in 
the plasma frame, which immediately dissipates its relative 
kinetic energy into heat. We will assume that this heat goes 
entirely to the proton component. The dissipated energy per 
decayed neutron is (F^e/ — l)mpC^, which gives the heating 
rate. 


dQp 

dt' 


(Frei - 1) nripC —. 


(69) 


Protons also exchange energy with electrons via 
Coulomb collisions with rate dQep/dt' (eq. 22). The net ther¬ 
mal balance of the proton component is then given by 


Our dynamical model of the neutron-loaded outflow may be 
summarized as follows. It is described by different equations 
before and after transparency, and the equations match at 
the transparency radius. Before transparency, r < Rt, we 
neglect the decay of neutrons. The evolution of Fp and F„ 
is found from equations lHHl and lIK^ . 

After transparency, r > Rt , we neglect the n-p collisions 
and assume F„ = const. We take into account the /3-decay 
which becomes increasingly important at larger radii. The 
evolution of Fp is described by equation 1651 . It includes the 
heating term dqp/dr which may be non-negligible at r ~ Rp. 
This term couples the dynamics with the proton thermal 
balance which is discussed in section 3.5 below. 

The two descriptions match at the beginning of the 
coasting phase = 0 where neither n-p collisions nor /3- 
decay affects the outflow. 


3.5 Thermal balance 


Thermal evolution of electrons and protons obeys the hrst 
law of thermodynamics leg. I16II . Using equations GU and 
we derive 


(1 + 0 


dr 


Ti 


(7-l)(l+0 


Ti_ _]_ d% 1 

n dr knp dt' c Fp /3p 


dqp dQ ep 

dt' dt' 

+ Fi.0i?T.7i (F,.e/ 1) nipC f ThptJnpC “h _ — 1(70) 

y J- relT/S J 

The /3-decay weakly affects the plasma thermal balance as 
long as (Fre.iTp)~^ <C UpOnpC. The decay becomes important 
at large radii, after the outflow becomes transparent. 

The thermal-balance equation for electrons is similar to 
the neutron-free case (section 2 ), 


dQe 3 (Tc — Te) dQep 


(71) 


3.6 Destruction of neutrons by inelastic collisions 

We assume in this paper that neutron richness ^ changes 
only as a result of /3-decay and neglect other channels of 
neutron conversion to protons. In fact, near the decoupling 
radius R„p the collisions between neutrons and protons be¬ 
come sufficiently energic to produce pions (DKK99a). Thus, 
some collisions are inelastic, which may destroy the neutron 
component. 

Inelastic n-p collisions may convert neutron to proton 
n + p^p + p + n~ as well as proton to neutron n -|- p —> 
n -I- n -I- tt’*' . These reactions have equal cross section which 
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has been measured down to the 140 MeV threshold (e.g. 
Daum et al. 2002). The reactions have equal rates and do 
not change neutron richness 

Inelastic n-n collisions may destroy the neutron compo¬ 
nent via reactions n-|-n^d-|- 7 r“ and n + n —> n + p + tt~ . 
These reactions have same cross sections as reactions p+p 
d + Tv^ and p+p —> n-|-p-|-7r''", respectively, which have been 
studied in experiments (Shimizu et al. 1982). The cross sec¬ 
tion of dn~ channel does not exceed 0.1 of the elastic cross 
section and is less important. The main reaction of neutron 
destruction is n-l-n —> n-|-p-|- 7 r“. Its cross section becomes 
comparable to the elastic cross section when neutrons collide 
with relative energy E„n > 700 MeV and quickly decreases 
at smaller Enn- 

The rate of elastic n-n collisions is ^ times higher 
than the rate of n-p collisions, which equals the expan¬ 
sion rate at decoupling. Therefore, the timescale of n-n col¬ 
lisions at r ~ Rnp is times shorter than the expan¬ 
sion time Texp, and at large ^ neutrons may be treated as 
Maxwellian gas. This gas is heated by n-p collisions which 
dissipate energy (Trei — l)mpC? 12 ~ per collision, 

and the heating rate of neutrons is equal to that of pro¬ 
tons lea. 1671 . The neutron temperature kTn = {2/3)En may 
be estimated by integrating the heating rate per neutron 
dEn/dr « {nipRl2r)Pxei « {nripR/2r){r/RnpY from r = 0 
to r = Rnp- This gives a modest value kT„ ~ 100 MeV. 
A more accurate calculation, which includes adiabatic cool¬ 
ing of neutrons, may give even lower r„. We conlcude that 
the mean relative energy of n-n collisions, Enn = SkTn, is 
well below 700 MeV, so only a tail of the quasi-Maxwellian 
distribution will contribute to destruction of neutrons. 

Therefore most of n-n collisions at r ~ Rnp are expected 
to be elastic, and a small fraction C, of these collisions will 
convert neutrons to protons. The lifetime of a neutron at 
r Rnp is Tn^p ~ C~^i~^Texp- An initially high neutron- 
to-proton ratio ^ will be reduced by conversion at decoupling 
if ^~^Tn^p < Texp, i.c. if the initial ^ exceeds 

The exact ( and the corresponding upper bound ^max 
may be found with detailed kinetic calculations of collisions. 
Such calculations are not attempted here, and we consider 
an optimistic range of ^ < 10 which might remain intact af¬ 
ter decoupling. The case = 4 is chosen as a main example 
model. 

3.7 Pair cascade initiated by 7r° 

Collisions between baryons near decoupling can produce 
neutral pions via reactions p + p—>^p + p + Tr^, n + p ^ 
d + 7v° and n + n^n + n + n°. Decay of 7 r° initiates a pair 
cascade in the outflow (DKK99a). 

The inelastic p-p collisions are relatively rare at ^ > 1, 
especially when protons are cooled by Coulomb collisions 
with electrons. Reaction n + p ^ d + has a maximum 
cross section ~ 0.1 of the elastic cross section, and yields 
only ~ 0.1 Tv° per proton at decoupling. 

The cross section of reaction n + n n + n + iv^ is ~ 1/4 
times smaller than cross section of n+n n+p+TT~ . (These 
reactions have the same cross sections as the experimentally 
studied reactions p + p p + p + TT^ and p + p —> p + n + iv'^, 
respectively.) Hence, the timescale of production by a 
neutron is 4 times longer than which is in turn longer 

than Texp (section 13.611 . This implies that up to ~ 0.2 7 r° 


is produced per neutron at decoupling, which translates to 
~ 1 tt'’ per electron. 

A typical produced has energy of a few hundred 
MeV. It immediately decays into two photons, and the 
high-energy photons are absorbed by the thermal radiation 
(DKK99a). As a result, a relativistic pair is created which 
upscatters thermal radiation, and the upscattered photons 
create a new generation of pairs. A maximum possible pair 
yield of the initiated cascade, Y « 0.1, would be achieved 
if the cascade were ‘saturated’, i.e., if the upscattered 7 - 
rays always got absorbed by softer photons (Svensson 1987). 
Thus maximum 10 per cent of the energy could possibly 
be converted into pair rest mass, which corresponds to ~ 30 
per injected The cascade is not, however, saturated: 
after a few generations, the upscattered 7 -rays are able to 
escape the thermal radiation (DKK99a). Therefore, a more 
realistic pair yield is a few per cent. It corresponds to roughly 
10 injected per proton in the outflow. Most of these pairs 
annihilate soon after injection at r ~ Rnp- 

The inelastic n-n collisions also happen at r > Rnp 
with increasing timescale oc oc r^, and the pair creation 
remains significant until r ~ 4i?„p. The created pairs in¬ 
crease the opacity of the outflow and reduce the timescale 
of Coulumb interactions between protons and electrons. In 
numerical models presented below this effect is neglected. 

3.8 Numerical models 

Figures 31 a-nd show two numerical models of neutron- 
loaded outflows with initial neutron-to-proton ratio ^0 = 4 . 
The first model has 77 = 300 and the second model — 
7 = 3 X 10®. In both cases, 77 > 77* and the decoupling 
Fp > r„ takes place. It is especially signihcant in the high-77 
model. 

The evolution of Lorentz factors r„ and Fp has three 
stages: (i) Acceleration of both components which ends with 
different saturated values of Fn and Fp. Note that neutrons 
decouple earlier in the 77 = 3 x 10® model and therefore have 
a smaller r„ compared to the 77 = 300 model, (ii) Coasting 
stage r„ = const and Fp = const- (iii) Deceleration of the 
proton component by decaying neutrons. At about the same 
time the outflow begins to experience deceleration by an 
external medium. We do not study the external deceleration 
in the present paper and stop our calculations at r = 3 x 
10 ®® cm. 

The evolution of proton temperature Tp has an initial 
decline followed by two pronounced peaks. The hrst peak 
happens at the end of neutron acceleration. It is caused by 
the strong friction between the n and p components near the 
decoupling radius Rnp- The second peak accompanies the 
deceleration of protons at r ^ Rp — it is caused by absorp¬ 
tion of decayed neutrons by the proton outflow. The overall 
thermal evolution is markedly different from the neutron- 
free case and we describe it below in more detail. 

At the beginning of outflow expansion, the thermal evo¬ 
lution is similar to the neutron-free case: all components are 
thermally coupled at a common temperature. The temper¬ 
ature is controlled by radiation (which strongly dominates 
heat capacity of the outflow) and decreases according to adi¬ 
abatic law with 7 = 4/3. 

The behaviour changes when the Lorentz factor Fp « 
r„ reaches a value about 50. Then the proton temperature 
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Figure 4. Evolution of a neutron-loaded outflow with = 4, 
T] = 300 , L = 10 ®^ erg s“^, and Rq = 10 ^ cm. Upper panel: 
Lorentz factors of proton and neutron components (solid and dot¬ 
ted curves, respectively). Bottom panel: Temperatures of protons 
(solid curve), electrons (dashed curve) and radiation (dash-dotted 
curve). Vertical dashed lines indicate characteristic radii Rgb, Rt, 
Rep, Rej and R/s (see the text). The decoupling radius Rnp is 
close to Rsb- 
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Figure 5. Same as Fig. |4]but for t; = 3 X 10®. 
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Figure 6. Evolution of outflows with the same parameters as in 
Fig. 4 but with different = 1 (dashed curves) and = 10 
(solid curves). The model with = 4 from Fig. 4 is shown by 
the dash-dotted curves. Upper panel: Lorentz factors of proton 
and neutron components. Lower panel: proton temperature. 


decouples from Te and Tr and begins to grow. This is the 
result of the frictional heating and the quickly increasing 
relative velocity between n and p components. Only a frac¬ 
tion of the frictional heat is kept by baryons and a quasi¬ 
steady energy circulation is maintained in the accelerating 
outflow: radiation ^ relative kinetic energy of the n and 
p components ^ baryonic heat —> electrons —> radiation. 
The thermal balance of protons in this circulation controls 
Tp, and Tp quickly grows as the outflow approaches the n-p 
decoupling radius. 

After the decoupling radius, the n-p collisions become 
rare and the frictional heating extinguishes. On the other 
hand, the protons get thermally decoupled from electrons 
because of a long timescale of Coulomb energy exchange. 
The subsequent decrease of Tp is mainly caused by adiabatic 
cooling. 

Adiabatic cooling continues until the heating by (3- 
decay interferes the thermal evolution. The timescale of adi¬ 
abatic cooling in the proton frame is rjcTp which is pro¬ 
portional to r during the coasting stage. By contrast, the 
timescale of /3-decay is constant. Therefore the heating of 
protons by decayed neutrons wins the adiabatic cooling at 
some point, and Tp begins to grow as Tp oc r. This growth 
continues until the outflow approaches Rp where most of the 
neutrons decay. The exponential extinction of neutrons im¬ 
plies that heating extinguishes at r ~ Rp. Then the outflow 
again cools adiabatically. 

The two main effects of neutrons on the outflow — heat¬ 
ing and deceleration of the proton component — are espe¬ 
cially strong at high and high p. This is illustrated in 
Figures 1^ and |7| which show models with = 1,4,10 and 
p = 300,3 X 10®. At high and/or p the proton density 
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RADIUS [cm] 

Figure 7. Same as Fig. but for outflows with r] = 3 X 10^. 


of the outflow is low and the n-p decoupling occurs early, 
leading to a large relative Lorentz factor Frei between the n 
and p components. The large kinetic energy of the relative 
motion is then available for dissipation. 

We also note that outflows with very high rj and 
quickly become transparent and most of their energy is car¬ 
ried away by thermal radiation. This is illustrated by Fig¬ 
ure 1 ^ which shows three components of the outflow lumi¬ 
nosity: radiation, protons, and neutrons. The radiation lu¬ 
minosity is given by L~f = dTrr^c/JpFp (|aTr) for r < Rt 
and remains constant after The kinetic luminosities of 
protons and neutrons are Lp = A-k R cPpV'^pp c? and = 
Atx RcPnT^PnC?■ The sum of three contributions remains 
constant and equals the total luminosity L = 10®^ erg/s. 
One can see that the outflow with = 4 and rj — 300 emits 
about 30 per cent of its energy at the transparency radius, 
and the outflow with 77 = 3 x 10® — more than 90 per cent. 

All numerical models shown in Figures 4-8 assumed that 
the outflow starts to accelerate at radius i?o = 10^ cm. This 
is a reasonable assumption if the outflow forms near a stellar- 
mass compact object and expands either spherically or con¬ 
ically with a constant opening angle 9j > l/Fp. It is likely, 
however, that the free conical expansion is proceeded by a 
collimation stage (as envisioned by, e.g., collapsar model of 
MacFadyen & Woosley 1999). During the collimation stage 
6 j (r) decreases and the jet is hardly accelerated inside the 
progenitor star until it reaches the break-out radius Ro and 
starts free expansion with 6 j{r) — const and acceleration 
Fp oc r. Then our spherical/conical model applies with a 
large effective Rq. 

To see how the outflow dynamics is changed by col¬ 
limation we have calculated numerical models with Rq = 
10^, 10®, 10® cm. The results are shown in Figure 9. When 
i?o is large, a given Lorentz factor is reached at a smaller 
density of the outflow. As a result, neutrons decouple at a 



RADIUS [cm] 

Figure 8. Contributions to the energy outflow rate from radi¬ 
ation (solid curve), protons (dash-dotted curve) and neutrons 
(dashed curve). The sum of the three luminosities L-y + Lp + 
is constant and equal to L = 10®^ erg The upper panel cor¬ 
responds to the model shown in Fig. and the lower panel — to 
the model shown in Fig. |5] 

smaller Lorentz factor F^ and the outflow develops a larger 
Frei, leading to strong heating. This trend is observed in Fig¬ 
ure 9. Outflows with large Ro also have smaller Fp and most 
of their energy is carried away by thermal radiation released 
at the transparency radius. 

We also note that explosions with large Ro have larger 
decoupling radii Rnp and Rt , and a smaller decay radius Rp. 
Therefore the two peaks of frictional and /3-decay heating 
are closer to each other, and the stage of adiabatic cooling 
between them shortens. At even larger Ro ^ the two 
peaks would overlap. 


4 DISCUSSION 

In this paper we developed the theory of relativistic neutron- 
loaded outflows. The presence a neutron component signif¬ 
icantly affects the early dynamics of GRB explosions. In 
particular, the plasma temperature is increased by many or¬ 
ders of magnitude, and this heating can compete with other 
heating mechanisms such as internal shocks or dissipation 
of magnetic fields in the outflow. 

The effects of neutrons are pronounced in outflows with 
?7 > 77 * given by equation whose typical value is several 
hundred. Then neutrons and protons develop a substantial 
relative Lorentz factor which leads to strong heating and 
momentum exchange when neutrons decay. We also note 
that neutron-rich outflows with very high 77 3 > 77 * lose most 
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RADIUS [cm] 

Figure 9. Evolution of outflows with the same parameters as 
in Fig. 4 but with different Rq = 10® cm (dashed curves) and 
Ro = 10® cm (solid curves). The model with Rq = 10^ cm from 
Fig. 4 is shown by the dash-dotted curves. Upper panel: Lorentz 
factors of proton and neutron components. Lower panel: proton 
temperature. 


of their energy at the photosphere: their thermal radiation 
carries most of the energy at the moment of transparency. 
Therefore, the neutron effects may be especially strong in 
GRBs with a significant thermal component in the radia¬ 
tion spectrum (a number of such GRBs have been identified 
recently, see e.g. Ghirlanda 2003; Ryde 2004). 

The calculations in this paper focused on the simplest 
model of a uniform hydrodynamic outflow with weak mag¬ 
netic fields. We did not consider, for instance, internal shocks 
(e.g. Rees & Meszaros 1994) and possible pair creation by 
nonthermal y-rays generated in the outflow. Pair creation 
may extend the optically thick stage of expansion, and the 
trapped radiation may convert its energy more efficiently 
into bulk kinetic energy of the plasma. A large-scale mag¬ 
netic field may gradually collimate the outflow so that the 
conical geometry of expansion does not apply (Vlahakis, 
Peng & Konigl 2003). 

The /3-decay is likely to affect the development of inter¬ 
nal shocks in the outflow. The shocks are caused by a non- 
uniform profile of the Lorentz factor, and the drag effect of 
decayed neutrons tends to smoothen this profile. The fastest 
portions of the outflow are more effectively decelerated and 
the initial contrast of Lorentz factors may be substantially 
reduced already at r ~ 10^“^ cm. This effect constrains the 
dissipation efficiency of the shocks, which is sensitive to the 
contrast of Lorentz factors (see Fig. 3 in Beloborodov 2000). 
In addition, the high temperature of the outflow heated by 
/3-decay may prevent development of the shocks. We defer a 
detailed study of these effects to a future work. 

The impact of neutrons on the prompt burst and its 
afterglow provides a unique opportunity to link the observed 


emission with physical conditions in the central engine of 
the explosion. The very presence of neutrons is a signature 
of an extremely hot and dense engine. Observable effects of 
neutrons may shed light on the mechanism of GRB trigger. 
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